Classical wave experiments on chaotic scattering 



U. Kuhlt, H.-J. Stockmannj, R. WeaverJ 

fFachbereich Physik der Philipps-Univorsitat Marburg, Renthof 5, D-35032 
Marburg, Germany 

^Department of Theoretical and Applied Mechanics, University of Illinois, 104 S. 
Wright Street, Urbana, Illinois 61801, USA 

Abstract. We review recent research on the transport properties of classical 
waves through chaotic systems with special emphasis on microwaves and sound 
waves. Inasmuch as these experiments use antennas or transducers to couple 
waves into or out of the systems, scattering theory has to be applied for a 
quantitative interpretation of the measurements. Most experiments concentrate 
on tests of predictions from random matrix theory and the random plane 
wave approximation. In all studied examples a quantitative agreement between 
experiment and theory is achieved. To this end it is necessary, however, to 
take absorption and imperfect coupling into account, concepts that were ignored 
in most previous theoretical investigations. Classical phase space signatures of 
scattering are being examined in a small number of experiments. 



PACS numbers: 03.50.De,03.65.Nk,43.35.+d 

Submitted to: J. Phys. A: Math. Gen. 

E-mail: ulr ich . kuhlOphysik . unl-marburg . de 
E-mail: r-weaver9uiuc.edu 



Classical wave experiments on chaotic scattering 



2 



1. Introduction 

Experiments with classical waves, originally conceived to study the spectral properties 
of closed systems, have more recently become an important tool to study the scattering 
properties of open systems as well. 

At about 1980 the main experimental information on spectra and scattering 
properties of chaotic systems came from nuclear physics and resulted in the 
development of random matrix theory by Wigner, Dyson, Mehta and others on one 
side and of scattering theory by Weidenmiiller and coworkers on the other side j^. 

Interest in the properties of chaotic systems suddenly renewed, when it was 
conjectured by Casati et al [J! and Bohigas et al 0] that the universal properties of all 
chaotic systems should be described by random matrix theory. An equivalent of the 
conjecture for open systems was formulated somewhat later by Bliimel and Smilansky 
0. These conjectures proved to be extremely useful and gave a new impetus to the 
whole field. 

The experimental situation, however, remained unsatisfactory until 1990, when 
the first analogue experiments were published, starting with vibrating solids |H] and 
chaotic microwave cavities 0. Meanwhile the technique has been extended to water 
surface and pressure waves |H1 E] and optical systems ^3 El ■ At about the same 
time the tremendous progress in the fabrication of semiconductor microstructures 
allowed to study the transport properties of submicron size structures of various types 
such as quantum dots tunnelling barriers and quantum corrals '15^. The 
state of the art till 1999 can be found in chapter 2 of reference |16j . 

In the beginning these experiments focussed on the spectral properties of closed 
systems though there was already an early experiment looking for the scattering 
properties of an open microwave system |17j . Strictly speaking, the systems are 
always open, due to the necessity to introduce antennas or transducers to perform 
the measurements. Therefore scattering theory is mandatory for a quantitative 
interpretation of the measurements. This is why it became necessary to establish 
a relation between the Green's function of the closed system and the scattering 
matrix |18| . which is a direct equivalent of a corresponding expression derived in 
nuclear physics many years ago. For the case of an isolated resonance the well- 
known Breit- Wigner formula is recovered jl9j . Details will be given in section 2. 
The same equivalence also holds for quantum-dot systems, as long as electron-electron 
interactions can be neglected. Thus with respect to scattering theory it does not make 
any difference whether atomic nuclei, quantum dots, microwave cavities, or vibrating 
blocks are involved, provided one is interested in the universal features only. 

This allows experimental tests of theoretical predictions from scattering theory 
with an unprecedented precision. A number of important consequences of scattering 
theory have been demonstrated with classical waves, but never found in nuclei, such 
as the algebraic decay of the scattering matrix from a system with a small number 
of channels |20l 12 1| . or the fact that with increasing coupling to the channels the 
resonances eventually become narrow again 22 , known as resonance trapping. 

There are a number of clear advantages of microwave resonators and vibrating 
blocks as compared e. g. to nuclei, but also to quantum dots. First of all, the geometry 
is perfectly known, and the coupling to the external channels can be calculated 
explicitly from the antenna geometry. The systems are clean, and if there are 
impurities, they are introduced intentionally. The geometry can easily be varied thus 
allowing level dynamic measurements of different types, which is impossible for nuclei. 
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and can be achieved only in a very restricted way in quantum dots. Further advantages 
arise from the fact that the relevant wavelengths in microwave resonators or vibrating 
blocks are of the order of mm to cm resulting in very convenient system sizes. Last but 
not least the experiments can be performed at room temperature since there are no 
dephasing processes as in quantum dots which have to be frozen out. This advantage 
is lost in microwave experiments using superconducting cavities, but still the applied 
temperatures of some K are far above the 0.01 to 1 K temperatures usually applied 
in quantum dot experiments. 

In this review we present results from microwave resonators, vibrating blocks, 
and optical systems, and their comparison with the predictions of scattering theory. 
Wave experiments on localization in disordered systems are excluded since they are 
treated in a separate contribution to this volume j28| . Results on quantum dots, 
tunneling diodes and other mesoscopic structures are excluded as well. They would 
have deserved a contribution of their own, and will be mentioned only in cases where 
a comparison with results from experiments with classical waves seems appropriate. 

This review is organized as follows. In the next section a Green's function 
technique is applied to establish a relation between the scattering matrix and the 
Green's function of the system. In addition the random superposition of plane waves 
model is introduced to describe amplitude distributions in chaotic cavities. Though 
demonstrated for the example of a flat microwave cavity where there is a perfect 
analogy to quantum mechanics, the results of this section can be applied without 
changes to all experiments with classical waves. Therefore this section is the basis for 
all following ones where the results in open microwave cavities, vibrating blocks, and 
optical systems are presented. 

2. The scattering matrix for billiard systems 
2.1. The billiard Breit-Wigner formula 

A microwave billiard is a special example for a scattering system, where the measuring 
antennas take the role of the scattering channels. We are now going to discuss the 
consequences of the external coupling for resonance positions and widths. 

Let us start with a two-dimensional quantum billiard of arbitrary shape, and 
assume that there are a number of channels coupled to the billiard with diameters 
small compared to the wave length. Let ip{r, k) be the amplitude of the field within 
the billiard while matter waves enter and exit through the different channels. '4'{r, k) 
obeys the Helmholtz equation 



Close to the channels the field is given by a superposition of incoming and outgoing 
circular waves. 



where rj is the position of the ith. channel, and hi are the amplitudes of the waves 
entering and leaving the billiard through the channel, and H^p [x) and h'^\x) are 
Hankel functions. Denoting the amplitude vectors of the waves entering or leaving the 
billiard by a = (ai, 02, . . . , cat) and b = (61, 62, ■ ■ • , ^Af), respectively, the scattering 
matrix for the billiard system is defined by the relation 




(1) 



ip{r,k) 



a^H^ihlr - nl) - b.H^ihlr ~ n]) , 



(2) 



b = Sa . 



(3) 
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The Oi and bi are obtained by matching fields and their normal derivatives at the 
coupling positions of the channels. This is achieved by means of Green's function 
techniques. Details can be found, e.g., in Section 6.1.2 of Ref . |^. Thus one obtains 

1 - iW^GW 

where the matrix elements Wnm of W contain the information of the coupling of the 
nth eigenfunction to the mth channel, and G is given by 

where 7J = — A is the Hamilton operator for the undisturbed system. 

An explicit calculation of the Wnm has to take into account the details of the 
coupling geometries. For antennas in microwave resonators this has been achieved in 
Ref. 121]. But it is common practice to treat the Wnm as free parameters which are 
taken from the experiment. 

An elementary transformation of equation yields the equivalent expression 

For the case of non-overlapping resonances and point-like coupling this reduces to 

Sij{k) = S^j - 2ijG{ri,rj,k) , (7) 
with the modified Green's function 

o(n.r„k)^j: (8) 

^ k^ ~ kl + ijY.k\i'n{rk)\ 
(see reference i^n are the real eigenfunctions of the closed systems and 7 describes 

the antenna coupling. The last equations constitute the billiard equivalent of the 
Breit-Wigner formula well-known from nuclear physics for many years 25 . Often the 
scattering matrix is studied on the time domain. To this end the Fourier transform is 
performed, 

S,,{t)^^ f S,,{k)e-'^Uu;, (9) 



where to = hk /2m. It follows 



Szj{t) = Sij S{t) - ^ X! '4'niri)ipnirj) exp 











i^lUJn H 





,(10) 



with r„ = {fij/m)J2k \i'n{rk)f- This holds for the quantum mechanical case. For 
non-dispersive classical waves we have instead lo = ck, and 7 has to be replaced by k^ 
in Eqs. Q and ||SJ). Thus one obtains 

Sijit) = 6{t) - 27c^ Vn(n)V'n(r'j)sin(w„i)e-^^"* , (11) 

n 

where now r„ = jcj^k l'0n('"fc)l^- In both cases F corresponds to the decay rate of 
the square of the field. 

Equation (|7|) shows that the complete modified Green's function G(ri,rj,k) 
can be obtained from a transmission measurement between two antennas of variable 
position |18j . A reflection measurement at one antenna only yields the spectrum 
and the modulus of the wave function |7||2S|- The 'true' wave function G{ri,rj, k), 
however, is never obtained since the resonances are broadened by r„. This is the price 
to be paid for every measurement. 
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2.2. The random plane wave approximation 

A model which proved to be very useful to describe field distributions and correlation 
in chaotic billiards is the random plane wave approximation. It assumes that the wave 
function can be described at any point, not too close to the boundary, by a random 
superposition of plane waves, entering from different directions, but with the same 
modulus of the wave number. 



V'(r) = ^ a(fc„)e"="^ |fc„|=fc. (12) 




The model has been independently introduced in different fields, among others in 
acoustics 27, and in quantum mechanics of chaotic systems 28 . It is known as Berry's 
conjecture in the quantum chaos community. The random plane wave approximation 
can be obtained from the Green's function of the billiard by averaging over a small 
energy window as was shown by Hortikar, Srednicki j29| and further exploited by 
Urbina, Richter |30j . 

As an immediate consequence of the central limit theorem one obtains a Gaussian 
distribution for the amplitudes of the wave function of closed systems, 

(13) 

where A is the billiard area (see e.g. Section 6.2 of Ref. ^16). In open systems real 
and imaginary part of ^p{r) are uncorrelated and both Gaussian distributed. 

The model allows the explicit calculation of various correlation functions, the 
most popular among them the spatial amplitude autocorrelation function 

c^(ri,r2) = {ip* {r:,)ilj{r:^)) 

^jMkr), (14) 

where r = — r^]. 

The random plane wave approximation can be used to calculate the distribution 
of line widths. The line width in a billiard with v attached channels is given by a sum 
over the squares of v Gaussian distributed random numbers, (compare the expressions 
following equations (|10|l and The distribution p^{x) is hence a distribution 



2 J T{v/2y 

The distribution of resonance depths for a reflection measurement can be obtained 
as well. According to equation © the resonance depths are proportional to |'0„(r)p, 
where r is the antenna position. The depth distribution function for a closed system 
is thus a distribution with v — \, 



..(x) = vAe-4^ (16) 

i. e. a Porter-Thomas distribution. For an completely open system one obtains a 
depth distribution with 1^ — 2^ which is just a single exponential, 

P2(x) = |e-4-. (17) 
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The model allows the calculation of averages of transients of various correlations 
of scattering matrix elements as well. From equation H1U|) we obtain for the average 

of 

C^At)^{\S^At)\')-[^ + ^) , (18) 

i. e. an algebraic decay is expected with an exponent reflecting the number of attached 
channels. For a very large number of channels the distribution concentrates in a 
narrow peak, which is again a consequence of the central limit theorem. In this limit 
there is a cross-over in the decay of Cij (t) to a single exponential behaviour, as is 
immediately evident from equation (|10|l . 

a.W^e-". (19) 

It follows that the Fourier transform of Cij{t), which due the the convolution theorem 
of Fourier theory is just the energy autocorrelation function of the respective scattering 
matrix elements, is described by a Lorentzian. This is the regime of the Ericson 
fluctuations well-known in nuclear physics for many years. It is noteworthy 

that the same universal fluctuation statistics were long ago noted in acoustics as well 

2.3. Acoustic and elastodynamic systems 

The expressions Q to (I19II of the previous subsections have direct analogs in acoustic 
and elastodynamic systems. The concepts of ray propagation, and the more abstract 
formalisms, are essentially identical. The differences are chiefly in the governing wave 
equation and in the details of experimental transduction. 

Like microwaves in 2-d billiards, waves in 2-d shallow water tanks |^ 03 are 
governed by equation with fc^ interpreted as uj'^/gh (g: acceleration of gravity, h: 
water depth). Chinnery W presents acoustic measurements in a deep water tank in 
quasi-2-d cylindrical geometries. These experiments differ from the corresponding 2-d 
microwave experiments in that they generally have non-Dirichlet boundary conditions 
at tank edges and lower Q's than do microwaves in billiards. A similar wave equation 
applies to the acoustics of a drum head, for which Teitsworth 35 has compared 
random matrix theory (RMT) predictions, for level statistics and wavefunctions. 
Capillary waves ^ have even lower Q's, and an interpretation of as {w"^ p/Ty/^ , 
where T is the surface tension. There is at least one report of wave chaos in a soap 
bubble |36| . for which the governing equation is again (0, with fc^ interpreted as 
uj'^T/ph, plus corrections for air loading. 

There are surprisingly few reports of wave chaos in 3-d acoustic systems. One 
example is seismic waves, i.e., low frequency elastic waves in the earth, are also 
governed by equation H2U|) . At short enough wavelength they are scattered by disorder 
and show mesoscopic behaviors. Further there are experiments in gasses, e.g., dry 
helium or argon, which are worth exploring, as impedance mismatches with boundaries 
can be made large and Q's correspondingly large also. It should be mentioned that 
there is a large literature on acoustics in rooms (e.g. IFt!, IHSf US] ) , however very little 
of it is carried out in a context of RMT. Notable exceptions are the works of Schroder 
EOlllIlllllllSl, Lyon 121], and Davy 

The chief acoustic venue for experimental studies of wave chaos and related 
mesoscopic wave physics has been elastodynamics, i.e., ultrasonics in solids. At MHz 
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frequencies wavelengths are of the order of several mm in typical solids, and laboratory 
systems are of convenient size. In aluminum, Q's can approach 10^, [ISIITTII^I^I^ . 
in quartz they approach 10^ |5U[ 1^ . The governing equation differs from ^ in that 
the dependent variable is a not a scalar, but rather vector displacement m, and A is 
replaced by a tensor-valued linear differential operator, £ij [S2| 



u)'^p{x)ui{x,u}) — £ijUj{x,uj) = ~p{x)fi{x,u)) (20) 
£i] = -dkCikji{x)di 

where / is the density of externally applied force divided by the mass density p, and 
summation convention is in force. In isotropic media £ is given in terms of the Lame 
moduli of elasticity A and /i, 

£ij = -{XSikSji + p5ij6M + pSiidjk)dkdi. (21) 

Propagation and reflection of rays differs from the simpler case of equation in 
that there are two propagation speeds corresponding the transverse and longitudinal 
polarizations (more in anisotropic media) and reflection generally results in mode 
conversion. A ray of one type incident on a wall will reflect into two or more types. 
Thus ray dynamics is complex. 

In 2-d, e.g., in a plate below the first cutoff frequency, and in the presence of 
up/down symmetry in material properties and in geometry the governing equations 
decouple into two classes of waves. One may be identified as flexural; it has a 
displacement that is predominantly normal to the surface. It has a governing equation 
- at sufficiently long wavelength - like Q with A replaced by A^, the Kirchhoff 
plate equation. These fiexural waves propagate with high dispersion. Except for the 
dispersion and the modified boundary conditions, such systems closely resemble the 
familiar scalar billiard. They have been discussed in the context of ray chaos and 
semiclassics and RMT by Bogolmolny and Hugues 53^. The other decoupled wave 
has displacement vector predominantly in-plane, and consists of superposed shear and 
dispersive-longitudinal motions, coupled by mode conversion at the boundaries, the 
Poisson plate equation. At frequencies above the first cutoff (/ ~ ct/2 x thickness) 
the elastodynamics of a plate is best described in terms of its several dispersive Lamb 
modes of propagation all coupled at the boundaries. 

Regardless of the precise form of the governing equation or the geometry, and 
indeed regardless of the presence of internal heterogeneities, all these systems may be 
abstracted in the limit of no absorption and in the presence of the proper boundary 
conditions, by a self-adjoint operator problem. With typical, e.g., traction-free, 
boundary conditions, the governing equation 1)20(1 remains not only Hermitian, but 
real and symmetric, i.e., time-reversal invariant, and the eigenvalues are real and 
non negative, and the eigenfunctions it" (a;) real and orthogonal. 



Ho\u"} ^ u;l\u"}; H, = {\/p)£ (22) 
= j u'\x)-u"'{x)p{x)dx = 5nm 

The mass density p acts as a weight function in the inner product. The Hamiltonian 
Hq is real-symmetric positive semi-definite. It is conjectured for many purposes to be 
ergodically equivalent to a member of the Gaussian Orthogonal Ensemble of random 
matrices. 
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Linear loss mechanisms are modeled by visco-elasticity, a general framework based 
upon a presumption that stress depends only on the history of local strain. In this 
case the otherwise Hermitian operator Hq gains an additional term, an imaginary 
symmetric operator —iWW'^, which may in principle depend on frequency. Such 
dependence is, however, expected to be weak. On expressing WW^ in terms of its 
(real) eigenvectors |c) and eigenvalues Xci one writes 



H = Hq ^ iWW^ ^ Ho - i^,\c)xc{c\ (23) 

This may be compared to the effective Hamiltonian in the denominator of 10). 
The eigenvectors c represent loss channels, the real quantities Xc their strengths. The 
channels |c) may be represented as normalised vector functions of position c(x). The 
Xc and |c) are expected to vary slowly with frequency. Some of the loss channels 
may be associated with physical contacts to the body, supports or attached sensors. 
Others may be internal and material. The non- Hermitian but symmetric operator H 
has complex orthogonal eigensolutions 

i7|ix") = (c^„-ir„/2)2|«") (24) 

It may be noted that the above inner product is defined without complex conjugation. 

Level widths are given exactly by (compare the expression following 
equation (fTT|l ). 

^ _ 1 ^ Xc|(t^"|c)P 

There is a similar identity for the natural frequencies 

2 Vl _ |(t."*|go|t.")| 

— - ^^^^ 
Green's function (cj^ — H)^^ is given by a modal expansion, cf equation (O, 



or 

i6"(a;)tt"(y) 
llP- - (u;„ - ir„/2) 



In the time domain one has, 

r<( T tt"(x)M"(y) exp{-ia;„t - r„V2} 

G X, y; t) ^ -Im > — — 29 

^ u;„-ir„/2 

<(x)«^"(y)exp{-iw„i-r„i/2} 
G,(.,y;t) = -Im^ ^„ - ir„/2 

Elastic waves are coupled into and out of the bodies by, most often, piezoelectric 
devices. The ideal device exerts a concentrated point impulse in response to an 
electrical impulse. In practice the dynamical force is distributed in time and over the 
face of the transducer. Similarly an ideal detector responds with an electrical signal 
proportional to the instantaneous displacement at a point. In practice a detector 
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gives a signal dependent on the recent history of displacement across the transducer 
face. Calculating the space-time transfer function is difficult. It is a complicated 
problem in piezoelectricity and elastodynamics |54|. It depends on how the transducer 
is coupled and on the material of the body. Nevertheless, general arguments using 
electro-mechanical reciprocity show that a detector d connected to a single coaxial 
cable may be represented as delivering a signal 

V = Xd{d\u) (30) 

linear in the displacement field. The detection channel (d| may be presented as a 
vector function of position, non zero only at the position of the physical detector and 
oriented in the direction of its sensitivity. Similarly a source s delivers a distributed 
force (normalised to unit mass density) f{x) in response to an applied voltage V, 
which may be expressed as 

l/> = \s)XsV (31) 

where |s) represents the spatial distribution and orientation of the force f{x) generated 
by the source V. \s) and \d) as well as the scalars Xd and Xs in general depend on 
frequency, but do so weakly on the scale of the level spacings and over the ranges of 
fluctuation expected in the wave response of the structures. 

Unlike the usual case with microwaves, our sources and detectors do not present 
significant compliance to the specimen. What impedance there is, is almost entirely 
mechanical and independent of the electronics. Signal measurements per se therefore 
do not absorb or affect the waves. While the contact devices may so affect the waves, 
it is usually found that the vectors \d) and \s) do not correspond to significant loss 
channels |c) . Thus unlike the microwave case we do not speak of measuring S matrices 
but rather of responses. The signal in a detector d due to an input F'" to a source s 
is the solution to 

H\u) - ,u'\u) = \f) ^ \s)xsV^- (32) 

= Xd{d\u) (33) 

i. 6., 

y-*=Xd(d|G|s)x.V- (34) 

The output signal is the input signal convolved with the ds component of the Green's 
function, convolved with the source and receiver functions Xs and Xd, respectively. It 
is for this reason that we speak of measuring components of the Green's function, not 
components of an S matrix. In acoustics literature the process is not usually termed 
scattering or measurement of an S matrix element, although it is scattering in the 
sense of this review article. 

The approximation of negligible loss into the channels d and s is exact if source 
and detector are non-contact. Laser interferometry provides a non contact, albeit less 
sensitive detection . Dropped ball bearings EHl 1^ i and laser thermoelastic 
excitation are also sources which do not act as loss channels. The use of contact 
pin-like piezoelectric transducers ^BJ 14 71 1551 15B| , or needles jSJI 021 is only minimally 
intrusive. 

The random plane wave approximation (|12|l does not apply directly to 
elastodynamics, as there is additional information needed on the relative amplitudes of 
the different wave types. Nevertheless, RMT tells us (for the non lossy case where the 
modes are real) that a diffuse wave field in general, and a mode in particular, should 
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be a Gaussian random process. The generalized Berry conjecture asserts that the 
two-point mode correlation function is the imaginary part of the Green's function 



where dN/duj is the spectral density of modes. 

This is equivalent to equation but applies also to non-scalar waves and to 
regions near heterogeneities. That this is exact if the average is over a frequency 
interval, i.e., for short times, was pointed out in [5BI . 

In the presence of losses, the modes are in general not real. But the above equation 
tells us that the real and imaginary parts are uncorrelated. That the imaginary parts 
are Gaussian processes also has been assumed |59j . but the relative variance of the 
real and imaginary parts is as yet not understood. 

As in the discussion preceding equation (|15l) and 125|) . one notes that a level 
width is given by a weighted sum of the squares of Gaussian random numbers. If all 
channels have equal strength x, and if the modes are approximately real (as they are 
in the weak absorption limit) then a line width F is a chi square random number, and 
|F°"*(t)P decays in time like an incoherent superposition of terms, each ~ exp(— r„t) 



where u is the number of channels and T is the average decay rate. This was first 
derived in acoustics by Schroder (60j and later by Burkhardt ES] . This differs 

from (fT^ due to the assumed lack, in acoustics, of significant losses into the source 
and detection channels. As in l|19(l . it reduces to a simple exponential in the limit of 
a large number of channels v. 

In section 0| we review the many experimental evidences in support of the 
applicability of RMT, for open (lossy) and closed (nonlossy) acoustic systems. Many 
of the earliest studies emphasized eigenstatistics, in particular level and eigenfunction 
statistics. These are reviewed for completeness. More recent research has emphasized 
responses, i.e., scattering. Response statistics are often measured with the intention 
to understand them in terms of eigenstatistics. Thus the two types of studies are 
closely related. 

3. Microwave experiments 

3.1. Fundamental tests of scattering theory 

There are a number of fundamental predictions from scattering theory on resonance 
shapes, width distributions, decay behaviour, and others, all of which have been 
verified in microwave systems. We give here only a short account, as this subject 
has already thoroughly been discussed in section 6 of reference |lfij . 

Alt et al jMJ performed a detailed study of the shapes of resonance in the non- 
overlapping regime and could show that they are described perfectly well by the billiard 
Breit-Wigner function (|7|) . In a quarter-stadium-shaped microwave billiard with three 
attached antennas the same group looked for the distribution of the line depths at one 
antenna, and the sum of the line depths at two and three antennas, respectively, and 
found distributions with v = 1,2,3 as expected In the time domain one 

expects for the system an algebraic decay with a power of 7/2 due to the presence of 



-lmGij{x,y ,ujn) 



dN/dio 
2uj 



(35) 




(36) 
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three antennas, see equation ((TH|l . which was also observed in [5n|. Algebraic decay 
behaviour has been observed in reverberant blocks as well (see Section ^3)). 

The first microwave scattering experiment has been performed by Doron et al 
| 17) in an elbow-shaped resonator. It will be discussed in more detail in section 13.61 
in the proper context. It is interesting to note, that essentially the same system 
was studied numerically independently by Jalabert et al jH7j but at this time in 
the context of universal conductance fluctuations, a nice illustration of the fact that 
nuclei, microwave billiards, and mesoscopic systems do not differ, as long as universal 
scattering properties are concerned. 

A measurement of the elements of the scattering matrix yield directly the 
billiard Green's function, though modified due to the presence of the antennas (see 
equation Q)- Thus the complete quantum mechanical information is available. This 
was used in reference |18| in a measurement in a quarter stadium billiard, where one 
antenna was fixed, whereas the other one was scanned through the billiard. From 
the transmission spectra the propagator is obtained by means of Fourier transforms. 
For short times a circular wave was found, emitted from the fixed antenna, being 
destroyed after a small number of reflections. After some time the primary pulse 
was partially reconstructed due to the focussing properties of the circular part of 
the boundary. Similar measurements have been performed in a ray-splitting billiard, 
realised by a rectangular billiard, where a quarter-circular insert of teflon was attached 
to one corner [HHj. Teflon has an index of refraction of n = 1.44, thus giving rise to a 
ray splitting upon every reflection. In microwave billiards the matrix elements Wnm of 
W, describing the couphng of the antennas to the wave function, see equation Q), can 
be explicitly calculated from the geometry [^21 1^ ■ Since the poles of the scattering 
matrix S are given by the eigenvalues of the effective Hamiltonian _ffcff = H — iWW'^ , 
see equation ©, this allows a detailed study of these poles in dependence of the 
coupling strengths Details will be given in section 

The distribution found for the line widths presents an indirect verification 
of the random plane wave approximation. A nice direct verification, though in the 
regime of visible light, has been obtained by Doya et al JOI- The authors studied 
the transport of a He-Ne laser beam through a glass fiber with a i?-shaped cross- 
section with R = 63.5 fim. For the near-field intensity directly at the exit of the 
glass fiber they found a speckle-like pattern. Introducing an additional lens, and 
mapping the intensity pattern in the focal plane, the far-field intensity was obtained. 
It corresponds to the modulus squared of the two-dimensional Fourier transform of 
the near-field amplitude. The far-field intensity was found to concentrate on a circle 
with no specific preference for any direction, thus demonstrating the validity of the 
random plane wave approximation. 

3.2. Pole distributions in the complex plane 

In the case of resonance overlap the Breit-Wigner formula is not valid anymore (see 
equation ©), and one has to consider the full S'-matrix as given in equation (|3J| 
resulting in new and partly counterintuitive behaviour such as the resonance trapping 
effect. 

Different scales of lifetimes appear for the resonant states if there is sufficient 
overlap j7()l 1711 17^ . Some of the states align with the decay channels and become 
short lived, while the remaining ones decouple from the continuum and become long 
lived, a phenomenon called trapping. The doorway states in nuclear physics provide 
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Figure 1. Resonance position for four resonances in a Sinai billiard with attached 
waveguide with variable opening (a) and their corresponding width Tn/'i versus 
the slit width d (b). The different resonances are marked by different symbols 
(taken from reference |22p. 



an example for the alignment of the short-lived states with the channels |73[ 17^ . For 
a clear experimental demonstration of the trapping effect, the coupling strength to 
the decay channels has to be tuned. This was done experimentally by Persson et 
al 1221, ^ detailed theoretical discussion can be found in [221 ■ In the experiment the 
microwaves were coupled to a resonator of the shape of a Sinai billiard via an attached 
waveguide, and the reflection was measured. Using a slit with variable opening width 
d the coupling from the waveguide to the billiard could be changed. The slit width 
was varied in steps of 0.1mm from 60 to 240 mm. Resonance positions and widths 
were extracted using a centered time-delay analysis In figure^a) the dependence 
of four resonance positions on the opening d of the slit is shown. In figure ^b) the 
width r is plotted as a function of the slit opening d. For small slit openings d the 
widths of the resonances first increase with d but for larger openings the width finally 
decrease again, thus demonstrating the effect of resonance trapping [221 • 

In the case of resonance trapping all resonances are treated on equal footings. If 
there is only an overlap of two resonances and all others resonances are far apart one 
can describe the occurring pattern in the transmission by Fano resonances [7S| [751 ■ 
Fano resonances have been observed originally in photoabsorption in atoms |77[l78ll7^ 
but occur in many fields in physics like electron and neutron scattering (HOI I81j . 
Raman scattering j82| . photoabsorption in quantum well structures |83j . scanning 
tunnel microscopy |84| . and ballistic transport through quantum dots (artificial atoms) 
[531 IHSl El IHHI ■ Close to a given resonance Ef the spectra can be approximated by 
the Fano form |7U1I5^ . 

lTiE,d)\^^\^^^M±^^mml (37) 

where d describes the coupling strength, E^{d) is the position of the ith resonance, 
Ti(d) its width, and qi{d) the complex Fano asymmetry parameter. Window 
resonances appear in the limit q while the Breit-Wigner limit is reached for 
l^l ^ 1 (see reference [SU] for details). It should be noted that, in general, q cannot 
be simply identified with the ratio of resonant to non-resonant coupling strength 
[ 891 191L 1^ . Moreover, since Fano resonances result from the interference between 
resonances related to the eigenmodes in the cavity, the parameter q depends very 
sensitively on the specific constellation of the involved resonance poles |7S1I75) . 
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Figure 2. Left and centre: Real and imaginary part of a spectrum of a half 
Sinai microwave billiard with absorption. On top the measured spectrum is shown, 
on bottom the reconstructed spectrum as obtained from the harmonic inversion. 
The resonance positions are marked by short lines. Right: Resonance widths for 
the displayed spectra in MHz (top), and normalised to the mean level spacing 
(bottom) (taken from reference IQSI t. 



In none of the previously mentioned experiments the Fano resonances could be 
systematically studied by varying the coupling strength. Again this has been achieved 
in a microwave billiard [HOI IHHl ■ The setup consisted of two commercially available 
waveguides which were attached centred both to the entrance and the exit side of a 
rectangular resonator. At the junctions to the cavity metallic apertures of different 
openings d were inserted. In the applied frequency range only two transverse modes 
for the rectangular cavity could be excited. A continuous transition from a narrow 
Breit-Wigner {q very large, d small) to a window resonance {q 0, d large) was found 
|90|. In a further investigation the effect of decoherence and dissipation on the Fano 
resonance lineshapes was discussed ^5] . 

Up to now only resonances with no or weak overlap or, as in the case of the 
Fano resonances, with just two strongly overlapping resonances have been investigated 
experimentally. For a complete exploration of the features of the scattering matrix an 
investigation of the poles in the strong overlapping regime is needed. By means of the 
harmonic inversion technique |M1 IHZI this now has become possible pS] . 

In figure l^^left and centre) the complex reflection coefficient of a half Sinai 
billiard is shown for the frequency range from 8 to 9 GHz (for details see |SH1)- It 
is compared with a reconstructed spectrum obtained from the harmonic inversion. 
The resonance positions are marked by small lines. Apart from a secular overall 
phase shift the reconstruction is nearly perfect. Global phase shifts are not serious 
as long as one is only interested in resonance positions, width and absolute values of 
the amplitude, and their behaviour in dependence of an external parameter. On the 
right (top) the widths are plotted versus the frequency and (bottom) versus energy, 
normalised to a mean level spacings of one. In the example the resonance widths are 
of the order of the mean level spacing A, but it was possible to penetrate into regimes 
where the resonance width amounted up to about ten times the mean level spacings. 
For isolated resonances one expects a Porter- Thomas distribution [HHl Q , whereas for 
strong coupling the tails of the distribution decay according to i'r(r) oc F"^ jlOOIIlOlj . 
This could be verified in the experiment j98| . 
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Figure 3. The distribution of the reflection Pfl(iJ) is shown for four different 
coupling strength (ra=0.116, 0.754, 0.989, 0.998) and total absorption (7=0.56, 
2.42, 8.40, 48). Additionally the theoretical curve (yellow line) and numerical 
results are plotted (blue crosses) (for details see I102II10SI V 




3. 3. Distribution of reflection and transmission coefficients 

The first experimental investigation of scattering tlieory was done by Doron et al 
who studied the reflection properties of an elbow system It concentrated on 

the influence of the classical escape behaviour on the reflection, and will be treated 
in subsection 13.61 in the proper context. It took about ten more years until the 
influence of coupling and absorption was investigated experimentally by several groups 
in microwave biUiards in more detail [103 [10112111011103 UHSI ■ Barthelemy et al 
developed a model to describe the coupling of a microwave antenna to the system 
and could separate their contribution Tc to the total resonance width T from the part 
To/im resulting from ohmic losses in the wall 24 . 

Mendez-Sanchez et al investigated reflection spectra in dependence of the coupling 
strength and the absorption |102| . The experiment was performed in a half Sinai 
billiard, where an additional absorbing wall could be introduced to increase the total 
absorption. They could verify the theoretical predictions for the limiting case of strong 
absorption |l()7j . The known result of weak absorption |l()8j could not be verified as it 
assumed perfect coupling which was not realised in the experiment. An interpolating 
formula for the distribution Pr^ of reflection coefficients in the case of perfect coupling 
was derived [TOS HOO] 

PrARo) = C,jf^^^ [Aa^/^-^ + B{1 Rf/^] , (38) 

with a = 7d/3/2, ^ = a (e" - 1), B = (1 + a - e"), Cp = [Ar{l + P/2,a)/a^ + 
Be~°' /a]~^ and the universality index is a parameter describing the dissipation 

P03 and r{x,a) = t^~^e~'^dt is the upper incomplete Gamma function. 
Meanwhile there exists an exact but lengthy expression by Savin, Sommers, and 
Fyodorov |11()| . The deviations from equation H38|l are typically small. From the 
distribution Pr,o the joint distribution of the reflection coefficient and the phase of 
the reflection coefficient {S — y^e'^) for incomplete coupling can be calculated via 

/I — (S)'^ \ ^ 1 
PnMS) = [j^^jj^,) ^PnARo), (39) 

where (S) = VT~^ and Rq = \So\^ with So = {S - {S))/{1 - {S)S). In figure[Slthe 
distribution of the reflection P{R) is shown for four different coupling strength Ta and 
absorption 7^. The theoretical curves shown in blue are obtained via a projection to 
the corresponding axis of PR,e{S). A good agreement between experiment and theory 
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Figure 4. (a) Distribution of the imaginary part of cavity impedance for two 
different couplings from 9 to 9.6 GHz. (b) Rescaled plot of two curves in (a) 
after dividing the experimental impedances by the impedances Zq of the antenna 
(taken from reference ll05l ). 



is found for P{R) |102[ 1103] . In a further work the phase distribution P(0), called 
the Poisson kernel, and the joint distribution were determined experimentally, again 
in perfect agreement with theory |l()Hj . For phase distributions see section WM 

Hemmady et al applied a different approach |1(J5[ ll(J6l lllll I112| . They did not 
study the scattering matrix but looked for the impedance matrix Z instead. These 
quantities are related via 



Zq describes the coupling port and is fixed by an independent measurement. The 
impedance matrix relates the applied voltages V to the currents at the antennas 
{V = ZI), whereas the scattering matrix relates the incoming flux with the outgoing 
flux. Using the impedance matrix Hemmady et al could remove the effects of coupling 
thus obtaining a normalised cavity impedance, which shows a universal behaviour 
|1()5I llllj . This is demonstrated by figure ^ 

Microwave graphs are another possibility to study the infiuence of absorption 
on the distribution of reflection coefficients, real and imaginary part of impedances 
etc. This has been done by Hul et al II 1 31 . Again a perfect agreement with RMT 
predictions was found. 

Though there are no analytical results available for the joint distribution of 
scattering matrix elements for more than one channel, apart from the distribution 
of the modulus of the eigenvalues jll4j . there are results for the distribution of the 
total transmission T. This quantity is of practical interest, since due to the Landauer 
formula the total transmission through an open quantum dot is, up to an universal 
vector, nothing but the conductance G, 



where the index n sums over all incoming and the index m over all outgoing channels. 
There are a number of predictions on the distribution p{T) of the transmission 
in dependence of the channel number for systems with and without time-reversal 
symmetry |115l 1116) . In the calculation the matrix elements had been assumed to be 
Gaussian distributed, and uncorrelated, apart from the constraint of the universality 



Z = Zo{l + S)/{l~S). 



(40) 




(41) 
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of the S matrix. The assumption of Gaussian distributed scattering matrix elements 
seems to be in contrast to expression (0} for the scattering matrix. It can be shown, 
however, that for chaotic cavities and perfect coupUng both approaches are equivalent 
| 117| (see section IIA.4 of reference jll8| for a discussion). Experimentally obtained 
conductance distributions for quantum dots were in serious disagreement with the 
calculations, though the influence of the break of time reversal symmetry could be 
qualitatively reproduced 119 . A possible problem arises from the fact, that the 
geometry of the dots defined by the gate electrodes is ill-defined, and thus the channel 
number as well. 

In microwave cavities this problem does not exist. Experiments have been 
performed in a chaotic cavity with two attached waveguides both on the entrance and 
the exit side. Time reversal symmetry was broken by inserting ferrite cylinders. The 
phase-breaking mechanism of the ferrite can be qualitatively understood as follows: 
Applying a magnetic field, the electronic spins in the ferrite perform a precession 
with the consequence that microwaves reflected from a ferrite surface experience a 
phase shift whose sign depends on the direction of propagation. This effect has been 
used previously by So et al |12()j and Stoffregen et al |121| to study time-reversal 
symmetry breaking in closed microwave billiards. However, the ferrite introduces 
strong absorption, and the phase-breaking effect is present only in a narrow frequency 
window in the wing of the ferromagnetic resonance. Thus the calculation of the 
transmission distribution has to be modified by introducing additional absorbing 
channels. If this is done, both the break of time-reversal symmetry and the channel 
number dependence can be well reproduced by the calculations |122| . 

5.^. Correlations of S matrix elements 

The Fourier transform of the spectral autocorrelation function 



In chaotic systems K{t) shows a hole for small values of t, resulting from the spectral 
rigidity. Extensive work on the spectral hole has been done, among others in nuclear 
physics 123^ and in the analysis of spectra from microwave cavities .124.^ . The spectral 
hole is thus a convenient indicator to test chaotic behaviour. 

In open systems the spectral autocorrelation function has to be replaced by 
scattering matrix correlation functions 




(42) 



yields the spectral form factor 




(43) 



C[S:„ S,,]{E) = {S:,{E) S,,{E + E)) - {S:,{E)) {SUE)) , 



(44) 



or, alternatively, by their Fourier transforms 



c[s:„ s^it) = {s:,{t) s^t)) - s{t){s:,{E)) {sum) , 



(45) 



where 




(46) 
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Figure 5. Normalized cross-correlation function C* [""tot ' ""tot! (^) ^ 

rectangular and (b) a chaotic billiard with C3 symmetry with the Fourier 
transform taken over the frequency range 13 to 14 GHz (the division by 
accounts for the antenna coupling, see the original paper for details). Dotted 
and dashed lines are theoretical results for GOE and POE, respectively, taking 
absorption into account. The dash-dotted line in (b) corresponds to the theoretical 
expectation for the C3 billiard (for details see the original work) (taken from 
reference I126I V 

The optical theorem 

a[:l=2il^ReSaa) (47) 

estabhshes a relation between the diagonal elements of the scattering matrix and the 
total cross section, often used in nuclear physics. Correlation functions of the diagonal 
elements of the scattering matrix may thus be alternatively expressed in terms of 
correlation functions of the total cross section |125| . 

As an example figure |S1 shows the cross correlation for a rectangular and the 
chaotic cavity, obtained from reflection measurements at two different antennas |126j . 
There was a strong absorption, which had to be considered in the theory. Whereas 
for integrable systems an exponential decay of the correlation function is expected, 
caused by the absorption, for chaotic systems the correlation is always close to zero, 
thus exhibiting the correlation whole. In both cases the experimental results follow 
the predictions. Discrepancies are found only for the rectangular billiard for small 
times. They are caused by the presence of the antennas making the systems pseudo- 
integrable and producing a small correlation hole which is absent in the ideal rectangle 

An experimental observation of the correlation hole in the regime of visible light 
has been achieved by Dingjan et al |11| . The authors studied a transmission spectrum 
of a He-Ne laser through an open resonator made up of two end mirrors, either curved 
or flat and a curved folding mirror. By misaligning the latter one very large aberrations 
could be introduced thus making the system chaotic. Taking the Fourier transform of 
the spectral autocorrelation functions both for the regular and the disturbed system, 
curves were obtained very similar to the ones shown in figure |S1 for the microwave 
system. 

Very recently the study of scattering matrix correlation functions has been 
extended to the subject of fidelity. The fidelity amplitude is defined as the overlap 
integral of some initial state with itself after the evolution under the infiuence of two 
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Figure 6. Fidelity amplitude for a chaotic microwave billiard. The smooth 
solid line shows the linear-response result, while the dashed line corresponds 
to the exact result. The perturbation parameter Acxp has been fitted to each 
experimental curve (see the original work for details) (taken from reference [12?! t. 



slightly different time evolution operators U{t), U'{t), 

f{t)^{m\uHt)u'{t)\m) (48) 

The modulus square of the fidelity amplitude, the fidelity F{t) = |/(t)P, had been 
proposed years ago by Peres to quantify the stability of a quantum-mechanical 
system against perturbations jl28j . The renewed interest in the topic results from 
its relevance for quantum-information processing. The fidelity amplitude as defined in 
equation (|48|l . is experimentally hardly accessible, as it requires the knowledge of the 
complete wave function. This was the motivation by Schafer et al |129| . to introduce 
a scattering fidelity, defined in terms of scattering matrix elements of the perturbed 
and unperturbed system. In the weak coupling limit the scattering fidelity approaches 
the ordinary fidelity H48|) . Perturbation theory predicts an exponential decay of the 
fidelity for t < tn, and a cross-over to a Gaussian decay for t > tn- tn — fi/ /S.E is 
the Heisenberg time, where /S.E is the mean level spacing |13(J| . 

In a microwave experiment the perturbation was achieved by moving one wall in 
a chaotic microwave resonator with two attached antennas. In addition averaging was 
performed by superimposing the results from different energy window and geometries. 
Figure El shows the result for different perturbations. For small perturbations the 
Gaussian decay dominates, whereas with increasing perturbations strength a cross- 
over to an exponential decay is observed. In the range accessible to the experiment a 
good agreement was found with a linear-response prediction of random matrix theory 
| 131| . Deviations of the linear response results from the exact results, obtained by 
super-symmetry techniques [1321 1133) amount to at most 10 percent in the regimes 
shown in the figure, which could not be detected in the experiment. 

Essentially the same results have been obtained in an experiment in elastodynamic 
billiards 11341 . They will be presented in section 



3.5. Current correlations 

A systems has to be opened if a measurement is to be performed. As a consequence the 
eigenvalues are turned into resonances and acquire widths. This has been investigated 
in the previous subsections. Now we want to investigate the effect of the openings 
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Figure 7. Intensity distribution Pp{I) for two different values of rigidity p. The 
insets show the corresponding intensity pattern I = for a quantum dot hke 
microwave billiard. 



on wave functions and currents. Due to the openness the wave function acquires an 
imaginary part, i. e., 

= i'r + iip,. (49) 

A convenient parameter to characterise the openness is the phase rigidity |pp of ip 
| 135| . where 

m + 

p—1 corresponds to the closed and p=0 to the completely open system. 

The random plane wave approach introduced in section 12.21 allows to calculate 
the intensity distribution. In open systems the amplitudes a{k) (see equation (|12|l ) 
obey the relation {a{k)a{—k)) = p{a{k)a*{k)). For a fixed value of p this implies a 
generalised Porter-Thomas distribution for the intensity, 
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(51) 



l-\p\\ 

In figure [7| the intensity distribution for a microwave billiard with attached leads is 
plotted. The transitional behaviour from a closed system to an open system is clearly 
visible [136| . Good agreement with the theoretical prediction is found, where the only 
free parameter p is taken directly from the data. 

For closed systems there is a large interest in nodal lines (see e. g. |13& ), whereas 
in open systems there exists only nodal points. The nodal points correspond to vortices 
of the quantum mechanical probability density 

j(r) ^Im[V'*(r)VV(r)]. (52) 

At the same time the vortices give rise to phase singularities which have been 
demonstrated in a rectangular microwave billiard with an introduced absorber ,135] . 
There are a number of theoretical results on the distribution of the modulus of the 
currents j or its components j^, jy | 14(J| . vortex spacing distributions |141l I142j . 
various spatial autocorrelation functions etc., using the plane wave approximation. 
In quasi-tow-dimensional systems there is a one-to-one correspondence between the 
quantum probability current and the Foynting vector '139' , allowing to determine the 
flow patterns experimentally. The left part of figure |S1 shows a typical example. This 
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Figure 8. Left: Typical probability current pattern. Centre: Experimental 
vortex pair correlation function g, together with the theoretical prediction. Right: 
vorticity spatial autocorrelation function (taken from reference [1571 ). 



has been used to verify these distributions and correlation functions by microwave 
experiments J-43, 137 . In the central figure [S] the vortex pair correlation function 
is shown, together with the prediction from the random plane wave approximation 
fnHITIT] . In the right part of the fi gure the spatial autocorrelation function of the 
vorticity w = (1/2) [V x ji] is shown, again together with the theoretical prediction. 
In both examples a good agreement between theory and experiment is found. 

Usually spatial correlations vanish for fcr — > oo. Therefore it came as a surprise 
when Brouwer showed that there exist a number of correlation functions with long- 
range correlations jl45| . i.e., with a surviving correlation in the limit kr —* oo. This 
is the case, e.g. for the connected correlator of the squared intensity, {I{r)'^I{r'Y)c 
showing a surviving correlation of varjpp for fcr — > cx) (the subscript "c" refers to 
the connected correlator, {AB)c = {AB) — {A){B)). The long-range correlation was 
experimentally verified for various correlators jl36) . 

3.6. Classical phase space signatures of scattering 

In the preceding section we discussed some universal properties of the scattering 
matrix. In the following a number of features shall be presented associated with 
individual system properties. The basis for the semi-classical treatment of scattering 
is an expression for the scattering matrix originally given by Miller |14()| . 

5(r,r',fc) =^zi;„e^^", (53) 

n 

where the sum is over all classical trajectories from r to r'. Wn is a stability factor, 
which can be calculated from the classical dynamics, and Sn is the classical action. 
For billiard system Sn is given by 5„ — hkln, where Z„ is the length of the trajectory. 
Equation H53|l can be looked upon as a special result of Gutzwiller's semi-classical 
quantum mechanics, expressing the Green's function in terms of classical trajectories 

EH- 

In billiards systems the Fourier transform of equation (|53|l , with k as the variable, 
gives the stability weighted length spectrum. This has be used by Kim et al jl48l I137| 
in a number of papers to test a conjecture that the transport in quantum dots is 
mainly promoted via scarred wave functions. In a soft-walled microwave billiard the 
same authors found evidence for dynamical tunneling, i. e. quantum-mechanical cross- 
talking between stable islands in phase space being classically separated |149j . 
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For a large number of open channels the scattering matrix in the time domain 
decays exponentially with a time constant r, resulting in a Lorentzian shape of 
the corresponding autocorrelation function with width F — 1/r, see equation (|19|) . 
Semiclassical theory identifies r with the classical escape rate [151)1 1151) . This was 
the issue of the first microwave experiment on scattering in an elbow-shape resonator, 
where, however, the number of open channels was not large but just one. Lu et al 
|152| studied the spectral autocorrelation function in an n-disc scattering system with 
n = 1,2,3, and found a perfect agreement between the observed width and the values 
expected from the classical escape weight. 

A characteristic feature of the n-disc system is a classical repeller associated with 
a bouncing ball orbit between two neighbouring discs. It follows from semiclassical 
theory that the decay of the scattering matrix is no longer single-exponential but given 
by a superposition. 



where 7„' and 7„ are real and imaginary part of the eigenvalues of the Frobenius- 
Perron operator describing the classical dynamics of the system. The corresponding 
autocorrelation function shows the so-called Ruelle-Pollicott resonances. 

These resonances have been studied by Pance et al |153j for the n-disk system 
for n = 1, 2, 3. In all cases the found resonances were in perfect agreement with those 
expected from the classical dynamics. 

Another scattering experiment showing clear fingerprints of the underlying 
classical dynamics has been performed by Dembowski et al |154j . The authors studied 
the transport of microwaves through a channel confined from above by a Gaussian 
and from below by an inverted parabola. The dominant structure of the associated 
Poincare section, taken at either boundary, is a stable island associated with the 
bouncing ball orbit along the symmetry line of the channel. This region is classically 
not accessible. The trajectories surround instead the stable island repeatedly in phase 
space, before they leave the channel either to the left or the right side. This revolution 
show up in a lighthouse-like oscillatory decay of the scattering matrix on the time 
domain. 

3. 7. Summary 

Up to about 1990 the study scattering was essentially a domain of theory. Noticeable 
exceptions were the spectra of compound nuclei giving rise to the development of 
random matrix theory 0. The situation changed completely with the emergence of 
different types of billiard experiments. After the first microwave study [7] there were 
numerous experiments with classical waves and quantum dots |16| . In microwave 
experiments, in contrast to quantum dots, there is a complete control over the 
geometric parameters. Therefore a quantitative agreement between experiment and 
random matrix predictions was achieved for all cases studied up to now. To this end it 
was necessary, however, to take absorption and imperfect coupling into account which 
was ignored in most previous investigations. Fortunately there has been quite a large 
number of theoretical studies in recent years to close this gap f |155| . this volume), 
motivated last but not least by the possibility to test the theoretical predictions 
experimentally. 




(54) 
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4. Experiments with sound waves 

4.I. Measurements of eigenstatistics in acoustics 

Manfred Schroder (40| appears to have been the first to note that RMT level statistics 
should apply to acoustic systems and could have implications for room acoustics and 
for structural vibration. Interestingly, in his thesis he used microwaves in order to 
conduct experiments that would be analogous to an acoustic system. Weaver ^1 ITT] 
was the first to measure RMT statistics in acoustics. He determined about 150 elasto- 
dynamic levels in each of three small aluminum blocks excited by impulsive forces and 
measured in the time domain with piezoelectric detectors. Acoustic systems often lend 
themselves naturally to time-domain measurements. Such measurement are related 
by Fourier transform to the frequency sweeps typically done on microwave billiards. 
EUegaard and several co-workers repeated these measurements, and extended them 
|48l 0n| ■ Here EUegaard et al measured ~ 300 resonances in each of several aluminum 
blocks with varying spherical octants of material removed from a corner. Sub- 
Poissonian level statistics were observed in the less irregular samples and ascribed 
to mode conversions. Statistics approached the GOE as more material was removed. 
The work was extended to quartz blocks [SU], of higher Q, with the identification of 
about 1400 levels. Again various sized octants of material were removed from the 
corner, thus allowing a study of the transition between a sub-Poissonian and a GOE 
case. Neicu and Kudrolli |156| also studied the evolution of eigenstatistics under shape 
deformations. 

Later work from the Copenhagen group emphasized studies of 2-d plate- 
structures, analogous to the microwave billiard experiments conducted by many others. 
At sufficiently low frequency or small thickness of the plate it has only flexural and in- 
plane modes, which may or may not be coupled depending on the breaking of up/down 
symmetry. Bertelsen et al |157| studied nearly 1000 levels of a Sinai-stadium plate 
and found good agreement with a Weyl-like estimate of level density, and, if up/down 
symmetry is broken, good agreement with RMT predictions for other level statistics. 
Andersen et al |4H| showed how variations in air-pressure on such plates allows to 
distinguish between flexural waves (with non-negligible losses to air) and in-plane 
modes (with low losses to air). Once distinguished in this manner, the level statistics 
of the two kinds of modes could be separately shown to be GOE. 

Level dynamics, under variation of a parameter, has attracted attention also. 
Not only can one study the evolution of statistics as a system changes its symmetry 
or integrability, RMT also makes predictions for the (scaled) distribution of level 
curvatures as a parameter is varied. Bertelsen et al |51| studied the 'Devil's spaghetti' 
formed by the level dynamics under variations in temperature. 

At a lower confidence level, due to lower Q systems and few modes identified, 
others also have confirmed GOE level statistics in vibrations and acoustics, [HSI QH] 

The spatial statistics of diffuse fields in general and modes in particular, have also 
attracted attention. There are a few measurements directly imaging mode shapes, and 
a few of more generic diffuse fields. More commonly such statistics are investigated 
implicitly by means of field correlations amongst a discrete set of points. 

4-2. Spatial correlations 

Ebeling and RoUwage et al |^ EZ] discussed measurements of diffuse waves in 
reverberation rooms and in 2-d ripple tanks and showed agreement with predictions 
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based on the random plane- wave superposition argument. In particular they showed 
that field-field correlations {ip{x)ip{x + r)) = Jo(A:|r|), and that intensity-intensity 
correlations were the expected {I{x)I{x + r)) — 2Jg (fc|r|) -|- 1. The particular virtue 
of ripple tanks is ready imaging of wavefields. Chinnery [5] has also presented pictures 
of (2-d) acoustic modes in water. Similar studies have been made on drum heads |35| . 
and capillary waves [H|. 

Schaadt et al |57| studied mode shapes and their statistics in a plate with flexural 
and in-plane modes. By scanning their needle- like piezoelectric sensor, they could map 
a component of the surface displacement field of any mode of their choice. By means 
of pressure variations, they could determine if the mode in question was flexural or in- 
plane (or a mixture if the plate lacked up-down symmetry.) They showed the expected 
1 + 2Jq intensity- intensity {/-/) statistics amongst the flexural waves, and an expected 
more complex relation for the /-/ statistics of the in-plane modes. The complexity is 
due to the mixture of wave speeds associated with such modes . 

Several works of late have reported theory and measurements of diffuse-field- 
diffuse field correlations with features beyond the well-known Jo behaviour. These 
features are somewhat weaker and more difficult to discern than the short-distance 
Jo, but they are non-universal and of some practical importance. While this work has 
generally been confined to correlations amongst generic fields, one imagines that they 
apply also to modes. 

Weaver USHllISni showed that the field-field correlation function in a reverberant 
body was the (imaginary part of the) Green's function. This is essentially equivalent to 
a correlation between S matrix elements {Sij{uj) 8^^(10)). The correlation reduces to Jq 
in the interior of a large 2-d single wave speed isotropic billiard, but has other features 
in more complex bodies, or near boundaries or scatterers. In particular it retains 
signatures associated with reflections from nearby heterogeneities. Weaver derived 
this from a modal perspective [1581 1159j using assumptions on modal statistics, but 
it has also been derived in other ways. Derode showed that this relation is 

closely related to time-reversed acoustics. An example, taken from [Ifilj and plotted 
in figure|3 shows how a conventional waveform, with sharp arrivals in the time domain 
related to strong reflections, emerges from correlations of a diffuse field. Campillo et al 
p.62, 1631 ll64| have used such correlations to map the earth's surface structure using 
seismic waves. The ability to retrieve the Green's function by passive measurements 
is gathering attention in ocean acoustics and seismology |162| I163| II 641 11651 . 

4-3. Fidelity 

In section the concept of scattering fidelity was introduced to account for the fact 
that fidelity as originally defined is not accessible to the experiment. It was shown 
that under the assumptions of ergodicity the two definitions of fidelity are the same 
|167l 1132] . Snieder introduced the notion of 'coda- wave intcrferometry' in seismology 
| 158) equivalent to scattering fidelity. Weaver and Lobkis independently studied the 
same idea in ultrasonics |134| where temperature changes were responsible for the 
changed Hamiltonian. De Rosny et al |169j looked at a similar effect with acoustic 
waves in a tank with moving scatterers. Snieder applied it to ultrasonics in rock 1^ 
and showed an effect from temperature variations and from the opening of a crack. 
Weaver and Lobkis |134| showed that the measurement could be used to evaluate the 
mixing rates of different waves. In particular they showed that a regular body had 
much lower fidelity than did a ray-chaotic body. They confined their quantitative 
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Figure 9. A transient diffuse field in a multiply scattering water tank (a) 
is autocorrelated (b). On averaging (b) over 118 distinct distant sources, a 
sharp arrival appears. The sharp feature appears at a time corresponding to 
a retroreflection from a nearby heterogeneity, as would be seen if the receiver 
were used in pulse-echo mode (taken from reference llbll 'l. 




Figure 10. The enhanced backscatter peak near the position of the source in 
a 2-d chaotic elastic wave billiard was imaged by de Rosny et al (taken from 
reference 
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Figure 11. Measurements of the free decay of intensity |S(t)p. It is found 
that the Breit-Wigner formula 1361 and its generalizations (dashed lines), often 
do not match measurements. A full RMT calculation (solid lines) does match. 
Differences are most pronounced when channels are well coupled (taken from 
reference 1211 '). 

evaluations to times short compared to the Heisenberg time. Later Gorin et al |170| 
showed that RMT theory for the time dependence on the scale of the Heisenberg time 
was consistent with Weaver and Lobkis's measurements over such times. 

Ribay et al have considered the effect of temperature variations on reverberant 
acoustic signals in water and in air and on flexural waves in plates (171j . To the extent 
that temperature changes in these systems merely change wave speed and volume, and 
do so homogeneously, the only effect of such changes is to dilate the waveforms. This 
trivial effect is not difhcult to account for. Because longitudinal and shear wave speeds 
have different dependencies on temperature, in elastic bodies there is also a nontrivial 
change of waveform |134| . 

Lu and Michaels jl72j have suggested that ultrasonic fidelity, after correction for 
temperature variations, may be useful as a nondestructive means for monitoring the 
condition of materials and built structures. Growth of a crack for example will degrade 
fidelity (168) . and do so in a manner depending on the size and position of the crack. 

Lobkis (unpublished) has measured the almost perfect fidelity of ultrasonic 
waveforms under the addition of a very weakly coupled loss channel, and shown 
extraordinary sensitivity to that addition, far greater than that exhibited by the energy 
dynamics (see below) 

4-4- Reflection and Transmission Distributions 

Coherent backscattering, i.e., enhancements of the diagonal elements of G or S, are 
most profitably studied in the time domain. Acoustics is thus a natural way to examine 
it. Coherent backscatter enhancements by factors varying in time from two and three 
are readily observed: 15^^^ ~ (2 to 3) x |S'y P {i ^ j). Lyon 03] was the first to predict 
enhanced backscatter in acoustics, followed much later by a prediction of its dynamics 
by Weaver and Burkhardt jlTHI using a modal perspective and assumed GOE modal 
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statistics. It was first observed in acoustics by de Rosny et al in a 2-d billiard (See 
fig. I10|l and by Lobkis and Weaver in a 3-d billiard |174j . Larose et al |175j observed 
enhanced backscatter with seismic waves in an open multiply scattering system, with 
potential for application to characterizing the earth. 

A system consisting of two or more weakly coupled cavities can exhibit an 
enhanced backscatter that is a form of Anderson localization; all elements Sij or 
Gij with i and j within the same substructure are enhanced, by arbitrarily large 
factors, relative to others. This was observed by Lobkis and Weaver |176| in a system 
consisting of two coupled reverberant elastic bodies. It has been analyzed theoretically 
by Weaver and Lobkis |176| and by Gronqvist and Guhr [r77| 

In a non-lossy system, responses are given simply in terms of real modes and real 
natural frequencies. Spectra 15*^^(0;) ^ consist of distinct resonances; interpretation of 
spectra is straightforward. In the presence of loss these resonances spread and at high 
enough loss, they overlap, thus complicating interpretation of power spectra. At very 
high loss where the modes are well overlapped, the statistics simplify; they become 
equivalent to Ericson fluctuations ; responses in the time domain can be understood 
as simply decaying coloured Gaussian random processes. Frequency domain response 
statistics are derived from that picture by Fourier transform. 

Statistics at intermediate overlap are not as well understood. Most acoustics 
has concentrated on the question of the variance of power transmission |S'ij(ct;)p for 
i ^ j, (as opposed to its distribution as in section IX^ These may be considered 
intensity-intensity statistics, and thus an extension of the field-field statistics of the 
previous section. Lyon |44j early appreciated that these random spectra were relevant 
to reverberation rooms and structural vibrations. He made theoretical estimates for 
the intensity variance by expressing responses in the modal form H28|l . and taking 
expectations of |Gp written in terms of sums of various products of four modal factors. 
The square of the intensity T = jGp transmitted from a source at y with polarization 
j to a receiver at x with polarization i is given, after a few minor approximations and 
simplifications, by the modal expansion |178j : 



Under a naive assumption of simple eigenstatistics and homogeneous damping, 
(real modes, fixed level width) and using expressions like H55|l. Lyon evaluated 
expectations of T and T"^. Waterhouse [123 and Davy [IHOl CHIl UBS discussed 
and extended the analysis of Lyon. Davy carried out extensive measurements in 
reverberation rooms in order to corroborate such predictions, and found some limited 
success. Weaver 183] pointed out that consideration of spectral rigidity neglected 
by Lyon and Davy would improve agreement, but that discrepancies would remain. 
Lobkis et al |178| re-developed the modal-expansion theory for variance of T using 
exact GOE level and mode statistics. They also incorporated the effect of a level 
width distribution and consequent decay curvature. For lack of any other model, they 
retained the assumption of real eigenmodes, with modal statistics identical to those 
of the GOE. 

Like that of Davy, their theory consistently overestimated variance relative to 
measurements on ultrasound in reverberant elastic bodies, see figure lT^ They ascribed 
the discrepancies to their neglect of the complexness of the eigenmodes. A detailed 
comparison with a direct supersymmetric RMT calculation for variance (e.g. |184| ') 
has not yet been made. 




ul{x)u-{y) {v^\x)uf{y)Y u\{x)u\{y) {u>l{x)u'^^{y)Y 
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Figure 12. The normalized variance of power transmission |Gijp for i ^ j, 
as measured (jagged line) in a reverberant aluminum block, and as predicted 
(dashed line) using a 'naive' theory like that of Lyon 1441 or Davy I180lll8nil82l 
but incorporating GOE statistics on eigenvalues and real eigenfunctions, and 
a distribution of level widths. It is found that the naive theory consistently 
overestimates the variance (taken from reference 11781 ). 

Langley et al |185l I186j have continued to study the problem of variance in 
statistical structural vibration, and have presented evidence for a consistent difference 
between the RMT value for inverse participation ratio (three) and that observed in 
their structures. They furthermore conjecture that this difference may explain the 
differences between measured and predicted power variances. 

4-5. Dissipation 

In dissipative systems, losses can be intrinsic to the material or the structure, or be 
due to explicit open channels. The first order effect of losses is to introduce a uniform 
exponential decay in the time domain, equivalently to add an imaginary constant to 
the frequency. It may be a nuisance in experiments, but the effects are trivial when 
viewed that way, and the theoretical questions it engenders are simple. All eigen- 
statistics are preserved. 

In practice losses are not perfectly diagonal and there is some interesting 
mesoscopic physics. Dissipation is often distributed nearly homogeneously, so 
deviations from the predictions of simple theory are small; they are nevertheless 
observable. Three major non-trivial topics have been addressed. The simplest is 
perhaps that of 'decay curvature' |187l 11881 IHTl IH^ (discussed below) in which a 
variance of level widths manifests as a non-exponential decay of intensity in the time 
domain as in eauation ll8l and 1361 In the frequency domain this manifests in field-field 
correlations (5'y(w)S'j* (w-t-fi)) (see section lO|l . There are also studies of transmission 
power variances (l^l^) — (|S'p)^ (see section and of complex modes (see section 
Oil . 

Decay curvature H3f)|l was derived in acoustics by Schroder and later by 
Burkhardt |61[l62ll^ : it is based on a Breit-Wigner expression for level widths H15|) . 

(It was discussed for microwaves in eqn. Burkhardt found good 

agreement between observed decay curvatures and the formula and furthermore 
found that v correlated with the number of sites of material damage in a reverberant 
body. Weaver and Lobkis 21 and Lobkis Weaver and Rozhkov jl78| also found that 
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the formula H36|l does a good phenomenological job of fitting observed decay curvatures 
over accessible dynamic ranges. This is may be surprising, as the assumption of equally 
strong loss channels is surely incorrect. It is furthermore unreasonable to suppose 
that real modal statistics assumed in the derivation of (|18|l and H36|) are correct. Thus 
Weaver and Lobkis |21| went on to show that, carefully scrutinized, decay curvature 
does not fit to the formula (jSEJ. They did find that a full supersymmetric RMT 
calculation does fit the observations (see figure [TT)l . 

It is possible to understand the difference between the Breit-Wigner picture of 
decay curvature, equation H36|l and H18|l . and the exact RMT result by noting that the 
actual modes of a lossy system are complex. Thus the conclusion that level widths are 
chi-square random numbers given by a sum of squares of equal-variance Gaussian 
numbers is incorrect. That sum must be replaced by a sum of v squares of the real 
parts of the modes and v squares of their imaginary parts, cf equation Thus 
the level width distribution ought be more narrow in the presence of complex modes, 
and the decay curvature ought be less pronounced. RMT indeed predicts less decay 
curvature than does Breit-Wigner, especially for very open channels; see for example 
figure 111! Microwave results on non-algebraic decay have already been presented in 
section rm 

4-6. Complex eigenmodes 

Non-real modes, and in particular their statistics, which have been discussed for 
microwaves in terms of rigidity (see section 13.4(1 are at present under-studied in 
acoustics as well. Lobkis and Weaver showed that modes in moderately damped 
structures (where they remain discernable as they are only moderately overlapped) 
can have significant complexity, and that their phases have a distribution in good 
accordance with theoretical guesses (see figure They found that observed 

complexities were sufiicient to account for differences like those of figure ^] between 
naive theory and measurements. If theoretical estimates may be obtained for complex 
modal statistics, perhaps in terms of level width variances jlt)4j . it may be that the 
modal-based estimates for power variances IH HHOI CHIl EHl HHSl EHl HHHl EEE] could 
be brought into agreement with measurements. The distribution of resonance phases 
shown in figure [Ob is presumably related to that of the Poisson kernel introduced for 
microwave systems in section l3.4l but the exact correspondence is not yet clear. 

4-7. Summary 

In acoustics the study of response statistics has developed, over decades, mostly within 
the context of eigenstatistics. Theoretical predictions for field-field correlations and 
for intensity correlations have been formulated in terms of assumptions about the 
modes, with some success. It is only very recently that RMT theory and better 
informed measurements have been applied and the various assumptions questioned. 
For response statistics, as for decay curvature, a direct RMT calculation, e.g., 
by supersymmetry can be expected to match experimental measurements more 
perfectly than will an eigenmode perspective. Nevertheless eigenmode perspectives 
[ 441 11801 11811 11821 15^ that make convenient albeit incorrect assumptions about 
eigenstatistics have had a good deal of success. That approach to decay curvature 
is in most cases sufficiently accurate, and is also adequate for field-field correlations. 
The eigenmode perspective for intensity variances has so far failed to accurately 
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Figure 13. A short section of the real and imaginary parts of a transmission 
through a reverberant lossy elastodynamic body (left). It is clear that these two 
modes differ by a substantial phase. A histogram of the phase of resonances in 
this body fits well to a model of complex Gaussian modes (right) (taken from 
reference 1591 ). 

match measurements. If, however, such perspectives were to be extended beyond that 
of the most recent :,178 , one imagines that they might achieve an accuracy approaching 
that of RMT. While a full RMT treatments of systems with loss, using for example 
supersymmetric |184j or Monte-Carlo methods, are in principle possible, it remains 
desirable also to have an approach in terms of the familiar and useful and simpler 
concepts of eigenmodes. The estimates of Lyon, Davy, Weaver, and Langley are close 
enough to success to encourage the hope that such an approach may become feasible, 
and that a modal picture can be developed to comprehend all aspects of response 
statistics. 

5. Summary 

It is not unusual in science that related concepts are developed independently in 
different fields without mutual knowledge of each other. This is particularly true for 
scattering theory. The example of Ericson fluctuations in nuclear physics on the one 
hand, and universal conductance fluctuations in mesoscopic physics on the other hand, 
and universal power transmission coefficient fluctuations in acoustics, have already 
been mentioned. Another example is the random plane wave approximation which 
was proposed independently in the quantum mechanics of chaotic systems j28| and in 
acoustics ^7\. Also the concept of fidelity has been independently developed in both 
fields, though using different notions. It has been the motivation for this review to 
illustrate the underlying concepts common to all types of classical waves, with special 
emphasis on microwaves and sound waves. 

R. W. was supported by the NFS and H.-J. St. and U. K. by the DFG. 
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